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The predictions of our previous analysis about possible low-mass (m x < 50 GeV) relic neutralinos 
are discussed in the light of some recent results from WIMP direct detection experiments. It is 
proved that these light neutralinos are quite compatible with the new annual-modulation data of 
the DAMA Collaboration; our theoretical predictions are also compared with the upper bounds of 
the CDMS and EDELWEISS Collaborations. 
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Searches for neutralinos at colliders have not yet 
reached the sensitivity required to place a direct lower 
bound on the neutralino mass m x . The commonly quoted 
and employed bound m x > 50 GeV is derived from the 
lower bound on the chargino mass determined at LEP2 
( m x ~ GeV) under the assumption that the U(l) 
and SU(2) gaugino masses Mi and M 2 satisfy the stan- 
dard relationship Mi ~ at the electroweak scale. 
This hypothesis is a consequence of the assumption that 
these mass parameters have a common value at the grand 
unification (GUT) scale. 

In supersymmetric models with R-parity conservation 
and no gaugino-unification assumption at the GUT scale, 
an absolute lower limit on m x cannot be derived from the 
lower bound on the chargino mass. Instead, it may be 
established by applying the upper bound on the Cold 
Dark Matter (CDM) content in the Universe, VLcdm = 
Pcdm/Pc in combination with constraints imposed on 
the Higgs and supersymmetric parameters by measure- 
ments at colliders and other precision experiments (muon 
2, BR(b — » s + 7)). This point was discussed in Refs. 
where a lower bound on the neutralino mass of 
about 6 GeV was established as a consequence of the re- 
cent 2(7 C.L. upper-limit QcDAih 2 < 0.131 obtained by 
the analysis of the WMAP data 0. 

In Refs. P, 0| we also analyzed the properties of light 
(m x < 50 GeV) relic neutralinos from the point of view 
of their detection rates in WIMP direct search experi- 
ments. Two main properties were pointed out: 1) di- 
rect detection rates for neutralinos of mass lighter than 



50 GeV reach levels within sensitivities of current WIMP 
direct detection experiments; 2) in the mass range 6 GeV 
< m x < 25 GeV the detection rates are predicted to fall 
in a very restricted range (this is at variance with what 
happens at higher m x , where the detection rates in an 
effective supersymmetric MSSM scheme (effMSSM) vary 
over decades Q). As shown in Refs. [H, Q the property 
(2) is a consequence the fact that the detection rate has 
a lower bound induced by the upper limit on flcDAih 2 . 

Recalling that, for neutralino-matter interactions, 
coherent effects systematically dominate over spin- 
dependent ones, the aforementioned properties (l)-(2) 
are conveniently displayed in terms of the quantity 
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scalar 



is the neutralino nucleoli 



scalar cross-section and £ is a rescaling factor between 
the neutralino local matter density p x and the total 
local dark matter density po: £ = p x /po- Follow- 
ing a standard assumption, £ may be taken as £ = 
min(l, Q x h 2 / {SlcDhih 2 ) mill ) • 

The supersymmetric model considered in the present 
paper is an effMSSM scheme at the electroweak 
scale, with the following independent parameters: 
AI2 , p, tan (3, ttia , iriq, mj, A and R = Mi/M 2 . Notations 
are as follows: tan/? the ratio of the two Higgs v.e.v.'s: 
tan/3 =< i?2 >/< Hi >, p is the Higgs mixing mass 
parameter, vtia the mass of the CP-odd neutral Higgs 
boson, Triq is a soft-mass common to all squarks, TOj 
is a soft -mass common to all sleptons, A is a common 
dimcnsionlcss trilinear parameter for the third family, 



A~ h = Ai = Arriq and Af = Am^ (the trilinear paramc- 
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ters for the other families being set equal to zero) . Since 
we are here interested in light neutralinos, we consider 
values of R lower than its standard value: i?GUT — 0.5; 
jito.ton#,%A^^;in^!ii)| ake R in the range: 0.01 - 0.5. In 



the scanning of the supersymmetric parameter space, 
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given by the analytic expression 



A. Bottino, F. Donate, N. Fornengo, S. Scope! (2003) 




10 20 20 40 50 60 70 80 90 100 
rrt (GeV) 



FIG. 1: Ratio of the HI content to the B content in the 
neutralino composition as a function of the neutralino mass 
m x . The upper branch of the solid line derives from neutralino 
diagonalization (see text), while the lower one is obtained £cr 



from Q x h 2 < (f2cDM/i 2 )max = 0.131. 



we use the following ranges of the MSSM parameters: 
1 < tan/3 < 50, 100 GeV < \n\,M 2 ,m^mj < 1000 GeV, 
sign(/x) = -1, 1, 90 GeV < m A < 1000 GeV, -3 < A < 
3. We impose the experimental constraints: accelera- 
tors data on supersymmetric and Higgs boson searches, 
measurements of the b — > s + 7 decay and of the muon 
anomalous magnetic moment = (g^ — 2)/2 (the range 
— 160 < Aa M • 10 11 < 680 is used here for the deviation of 
the current experimental world average from the theoret- 
ical evaluation within the Standard Model). We notice 
that to satisfy the b — ► s + 7 constraint is not trivial at 
small tua values, since it requires some cancellation be- 
tween the contributions due to the Higgs-quark and the 
chargino-squark loops. However, it turns out that the ex- 
tent of compensation between the two terms is limited to 
an effect of 30-60 %. The range used here for the branch- 
ing ratio is 2.18 x 10~ 4 < BR(b -> s + 7) < 4.28 x 10~ 4 . 
For (flcDMh 2 )min we use here the value (QcDMh 2 ) m i n = 
0.095 [2j, derived at the 2a C.L. from the analysis of Ref. 

I- 

As discussed in Refs. [1-2], neutralino configurations 
at small m x have a dominant bino component with a 
small mixture with H®, i.e. writing the neutralino as 
X = aiB + a 2 !V (3) + a 3 H^ + a 4 .ff|, one has |oi| >> 
03 >> |a.2 1 , |«4 1- In this regime the ratio |<X3 1 | is 



P4 ~ sm6 w sin/3 ^ < 0.42 sin/3, (1) 

where in the last step we have taken into account the 
experimental lower bound /j, > 100 GeV. The allowed 
range of the ratio 1 0.3 1 / 1 a. 1 1 for increasing values of m x is 
displayed in Fig. 1. The upper boundary line is given 
by intrinsic properties in the diagonalization of the neu- 
tralino mass matrix and is obtained by maximizing the 
parameters M2 and tan /3 and by minimizing [i (for very 
small m x , Eq. (fT|) applies). The lower boundary line is 
derived from the condition fl x < (QcDM)max (using Eq. 
(5) of Ref. [2]). These two boundary lines fit well with 
the scatter plot obtained by a numerical scanning of the 
supersymmetric parameter space, also displayed in Fig. 
[T] For m x > 20 GeV the cosmological bound is satis- 
fied by the f-exchange in the annihilation cross-section, 
and 1 03 1 is no longer constrained from below. The lower 
bound on |ci3 1 / 1 a-i | in the scatter plot of Fig. [T]is due to 
the upper value (1 TcV) in the range of fi employed in 
our calculation. 

In Fig. [5] we show a scatter plot of the quantity 



(nucleoli) 
scalar 



as a function of m x . This scatter plot shows 



that, in the mass range 6 GeV < m x < 25 GeV, the quan- 
tity C^caiar 00 "^ f aus m a narrow funnel (see property (2) 
above); this funnel is delimited from below by configura- 
tions with f2 x /i 2 ~ (OcDM^ 2 )max = 0.131, and delimited 
from above by supersymmetric configurations with a very 
light Higgs boson (close to its lower experimental bound 
of 90 GeV) and with an VL x h 2 below (ricr>M^ 2 )min- For 
m x < 10 GeV only values of 30 < tan/? < 50 and 100 
GeV < < 300 GeV contribute, while in the interval 10 
GeV < m x < 25 GeV tan f3 extends also to lower values 
around 8 and is not significantly constrained. More- 
over, for m x < 20 GeV, tua is strongly bounded from 
above by {QcDAih 2 ) max , as shown in Fig. 3 of Ref. 0]. 
Notice that the dip at ~ 45 GeV is due to the Z-pole in 
the annihilation cross-section. 

It is also remarkable that, within the funnel, the size of 
£ CT scaiai COn ^ i s l ar g e enough to make light relic neutralinos 
cxplorablc by WIMP direct experiments with the current 
sensitivities. To illustrate this point, let us turn now to 
a comparison of our predictions with experimental data 

5, 6j which became available after our analysis of Refs. 

1 , J and with those of Ref. 0] ■ 

In Refs. 0, 0] we anticipated that a detector with a 
high exposure and a low threshold such as DAMA Q 
might provide significant information, not only for neu- 
tralinos with m x > 50 GeV, but also for neutralinos in 
the mass range 6 GeV < m x < 50 GeV. In Ref. 0| 
we could only give an estimate of the expected effects 
in the measurement of the annual modulation variation 
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FIG. 2: Scatter plot of ^scalar versus m x . Crosses (red) and dots (blue) denote neutralino configurations with Q x h 2 > 
(Hci3Mft 2 )miii and f2 x /i 2 < (r2cDA//i 2 )min, respectively ((f2cDA/ft 2 )min = 0.095) (a) The curves delimit the DAMA region where 
the likelihood- function values are distant more than 4a from the null (absence of modulation) hypothesis [jj; this region is the 
union of the regions obtained by varying the WIMP DF over the set considered in Ref. 0| . (b) The solid and the dashed lines 
are the experimental upper bounds given by the CDMS @| and the EDELWEISS [?| Collaborations, respectively, under the 
hypothesis that the WIMP DF is given by an isothermal distribution with a standard set of astrophysical parameters. 



performed by the DAMA Collaboration, since no anal- 
ysis of the experimental data at low WIMP masses was 
available at that time. 

Now, the recent presentation of new results by the 
DAMA Collaboration Q allows us to compare directly 
our theoretical predictions to actual experimental data. 
In fact, one has now the results of a much larger exposure 
than in the past, about 108,000 kg • day and, most im- 
portant, an analysis of the full set of experimental data in 
terms of a spin-independent effect over an unconstrained 
range for the mass of a generic WIMP. The results of 
this analysis are reported in Fig. [DJa), where the con- 
tour line (after Fig. 28 of Ref. [5() delimits a region of 
the m x — £o~j^ilr P iane > where the likelihood-function 
values are distant more than 4tr from the null (absence 
of modulation) hypothesis. In deriving this contour line, 
the DAMA Collaboration has taken into account a rather 
large class of possible phase-space distribution functions 
(DF) for WIMPs in the galactic halo. The categories of 
DFs considered in Ref. Q are those analyzed in Ref. 

the annual-modulation region displayed in Fig. [UJa) 
is the union of the regions obtained by varying over the 
set of the DFs considered in Ref. j9j|. From Fig. [Ha) 
we derive that the entire population of relic neutralinos 
with m x < 25 GeV as well BjS z\ significant portion of 
those with a mass up to about 50 GeV are within the 



annual-modulation region of the DAMA Collaboration. 
Thus, this yearly effect could be due to relic neutrali- 
nos of light masses, in alternative to the other possibility 
which we already discussed in Refs. 0] on neutralinos 
with masses above 50 GeV, and which is reconfirmed by 
the present analysis. 

Another experiment of WIMP direct detection, run 
by the CDMS Collaboration, has recently published new 
data Their results are given in terms of an upper 
bound 1 on a^, }^ for a given DF (an isothermal dis- 
tribution) and for a single set of the astrophysical param- 
eters: po = 0.3 GeV • cm -3 , v = 220 km • s _1 (i> is the 
local rotational velocity). This upper bound is displayed 
in Fig. [DJb) together with our theoretical scatter plot. 
Thus, we see that in case of an isothermal DF with the 
representative values of parameters given above, a size- 
able subset of supersymmetric configurations in the mass 
range 10 GeV < m x < 20 GeV would be incompatible 
with the experimental upper bound (together with some 



We consider here only the upper bound derived by the CDMS 
Collaboration without neutron subtraction. The upper limit ob- 
tained by subtracting an estimated neutron background appears 
too model-dependent; this point will be overcome, when the 
CDMS experiment is run in an underground location, as fore- 
seen by the Collaboration. 
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configurations with m x > 80 GcV). However, this con- 
clusion cannot be drawn in general. In fact, to set a solid 
constraint on the theoretical predictions, it is necessary 
to derive from the experimental data the upper bounds 
on £ CT scaiai COn ' ) f° r a large variety of DFs and of the cor- 
responding astrophysical parameters (with their own un- 
certainties) ; the intersection of these bounds would pro- 
vide an absolute limit to be used to possibly exclude a 
subset of supcrsymmetric population. An investigation 
by the CDMS Collaboration along these lines would be 
very interesting. 

Among other experiments of WIMP direct detection, 
the EDELWEISS experiment has published an upper 
bound which somewhat approaches the region of the low- 
mass neutralino population. This upper limit, again pro- 
vided for a single DF (the isothermal sphere with a stan- 
dard set of astrophysical parameters) is also displayed 
in Fig. HJb); it turns out to be marginal for the low- 
mass population, since it is tangent to our supersymmet- 
ric scatter plot (at m x ~ 30 GeV). The argument given 
before applies again in this case; one should vary the 
analytical forms of the DF, in order to derive a model- 
independent bound on C^calax • As for the low-mass 
configurations (m x < 50 GeV), since the current upper 
limit is already marginal for the isothermal DF, one does 
not expect any model-independent constraint. However, 
useful constraints could be derived for higher masses. 

In conclusion, we have shown that the experimental 
exploration of the low-mass neutralino population, theo- 
retically analyzed in our papers of Refs. 0, 0], is already 
under way in case of some experiments of WIMP direct 
detection and within the reach of further investigation in 
the near future. 

We have compared our predictions with available re- 
sults of various experiments separately, since the experi- 
mental results of different Collaborations arc not all de- 
rived under the same assumptions on the WIMP phase- 
space distribution function. A more effective compari- 
son of theoretical results with experimental data will be 
feasible, only when the analysis of different experimen- 
tal results in terms of m x — £^icalar° i s presented for 



each analytic form of the DF, separately. This is also the 
unique way of comparing results of different experiments 
among themselves. 

We finally notice that, in direct detection experiments, 
lighter WIMPs have to be faster, as compared to the 
heavier ones, in order to deposit a recoil energy above 
the energy threshold. As a consequence, for light WIMPS 
the calculation of expected rates and the determination 
of upper limits on the cross-section are very sensitive 
to the value assigned to the escape velocity and, more 
generally, to the details of the high-velocity tail of the 
DF. This introduces an important uncertainty, since for 
high-velocity WIMPS the assumption of thcrmalization, 
which for instance is assumed in all the models considered 
in the analysis of Ref. Q , is less robust than for the bulk 
of the distribution: non thermal components, such as 
streams, could have a sizeable or even dominant weight, 
affecting the usual estimates for expected rates. 

Note Added. 

After submission of the present paper, updated evalu- 



ations of the e + 



-based and r-based lowest-order 



(LO) hadronic polarization contribution to the muon 
magnetic moment have been presented by M. Davier et 
al. (hep-ph/0308213). In order to obtain a conserva- 
tive range for the deviation of the standard model value 
of a M from its experimental determination, we consider 
the e + — e~-based LO hadronic contribution by K. Hagi- 



wara et al. (hep-ph/0209187) together with the r-based 
LO hadronic contribution by M. Davier et al. Com- 
bining the evaluation by K. Hagiwara et al. with the 
other standard model contributions, we find for the de- 
viation from the world average experimental result the 
2<7 range: 133 < Aa M < 569. Instead, using the r— 
based Davier et al. result, one obtains the 2a interval: 
— 142 < Aa M < 286. If we combine conservatively these 
two determinations we finally have —142 < Aa M < 569. 
Employing this range instead of the one used in the 
derivation of the scatter plot in Fig. [2] does not mod- 
ify the features of the plot in any significant way. 
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